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Abstract 

Using the Yale stellar evolution code, models of fi Her based on asteroseismic 
measurements are constructed. A \ 2 minimization is performed to approach 
the best modeling parameters which reproduce the observations within their 
errors. By combining all non-asteroseismic constraints with asteroseismic mea- 
surements, we find that the observational constraints favour a model with a 
mass of 1.00io!o2 M <3i an a S e * = 6.433 ± 0.04 Gyr, a mixing-length 
parameter a — 1.75 ± 0.25 , an initial hydrogen abundance Xi — 
0.605to'oo5 and metal abundance Zi = 0.0275+°'.ooi- A* Her is in post-main 
sequence phase of evolution. The modes of I = 1 show up the characteristics 
of avoided crossings, which may be applied to test the internal structure of this 
type stars. Asteroseismic measurements can be used as a complementary con- 
straint on the modeling parameters. The models with mass 1.00 - 1.10 M 
can reproduce the observational constraints. Existing observed data 
of [i Her do not rule out these models. 
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1. Introduction 



Solar-like oscillations have bee n confirmed for several mai n-sequence and 
subgiant stars, suc h as a Cen A dBecdine . 2004T ). a Ce n B (Kieldsen et all 



120051) . Procyon A (jEggenberger et all 



2004bl ). T) Bootis (jCarrier et all 1200 



etc. The large and small frequency separations of p- modes can prov ide a good es- 
timate of the mean density and age of the stars (|Ulrich I . [l986l [l988). It has been 
proven that asteroseismolog y is a powerful tool f or determining the fundamen- 



tal p arameters of the stars (jEggenberger et al 
20061) . 



200,4 lEggenberger k. Carrien 



H Her (HR 6623, HD 161797, HIP 86974) is a G5 IV subgiant star, in the 
neighbourhood of the Sun. It is considered to be a solar-type star with mass 
1.14 MfT, effective temper ature T e = 5596 ± 80 K and r adius R = 1.77 ± 
0.07 R Q (|Fuhrmannl . ll998h . Recently, feonanno et all (|2008h detected solar-like 
oscillations on /i Her and identified individual frequencies in the range of 900 - 
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1600 fiHz. These seismological data will provide a constraint on the fundamental 
parameters of /i Her. In this work, we try to determine modeling parameters 
of n Her based on asteroseismic constraints using the Yale Rotation Evolution 
Code (YREC7) in its non-rotating configuration. 

The observational constraints available for /i Her are summarized in Sect. 
2, while the details of the evolutionary models and the computational method 
are given in Sect. 3. The results are presented in Sect. 4 and the conclusion is 
given in Sect. 5. 



2. Observational constraints 



2.1. Effective temperature, luminosity and chemical composition 

The effective tempe rature of ji Her giv en bv iFuhrmann (1998) is 5596 K, 
howe ver that given by Ivanov et al.l ( 2004) is only 5390 K. Comb ining other 
data (|Valdes et all 120041 iTakeda et all 12061 ISoubiran et~all . l2008h . we adopt 
an average effective temperature T e ff = 5500 ± 90 K. 

The luminosity of a star can be obtained through combining the knowledge 
of the magnitude and distance. By combining the visual magnitude V = 3.4 17 ± 
0.014, the bolometric correction BC = - 0. 15 ± 0.05 mag ( Fuhr mannl . 119981 ) . the 
solar absolute magnitude M bo i tQ = 4.746 ( Leieune et all 1998f) a nd the newest 
Hipparcos parallax n = 120.33 ± 0.16 mas ( Van Leeuwenf 2007 ) . we obtained 
a luminosity for /i Her of L = 2.70 ± 0.16 Lq 



In the version of the Catalogue of [Fe/H] determinations given bv lCavrel et al 
(l200lh . there are five metallicity values for /i Her. R ecent determination s give 
the values: 0.26 (|Soubiran fc Girardl . Exjj ) and 0.29 jTakeda et all l2005h . We 
adopt the average of these determinatio ns, [Fe/H] = 0.21 ± 0.07. This value 
is close to the value of 0.23 given by ( Fuhrmannl . 19981 ). For Population I 
stars, the ratio of surface heavy elements to hydrogen abundance is related to 
the Fe/H by [Fe/H] ~ log(Z/X) s - log(Z/X) Q , where (Z/X) Q is the ratio 
of the solar mixture. The most recent ratio of the heavy-el ement abundanc e 



to hydrogen abundance of the Sun, (Z/X) & , is 0.0171 (jAsplund et al.l . 120041 ) 



There a re, however , some discrepancies between this new value and seismical 



results (|Yang fc Bi, 2007). Thus in this work, we adopt the old value, 0.0245 



(|Grevesse fc Noeld . ll993h . Consequently, the value of (Z/X) s for \i Her is about 
0.040 ± 0.006. 



2.2. Asteroseismic data 



Bonanno et al 



sol ar-like oscillations for this star have been detected by _ 
( 2008[ ) with the SARG echelle spectrograph. Twenty oscillation frequencies have 
been identified by using modified and standard extraction methods between 900 
and 1600 ^Hz. By means of a least square best fit with th e asym ptotic relation 
of frequencies for all the identified modes, iBonanno et al. (2008) gave that the 
most likely value of the mean large and small frequency separation is Ai/ = 56.50 
± 0.07 /iHz and 6v = 5.03 ± 0.94 /iHz, respectively. Some of the observational 
constraints for /i Her used in this work are given in Table [TJ 



2 



Table 1: Observational data for /i Her used in this work. 

n [mas] 120.33 ± 0.16 

V [mag] 3.417 ± 0.014 

L/Lq 2.70 ± 0.16 

T eff [K] 5500 ± 90 

[Fe/H] s 0.21 ± 0.07 

At* [/xHz] 56.50 ± 0.07 

5v [/xHz] 5.03 ± 0.94 



3. Stellar models 



3.1. Input physics 

A grid of stellar evoluti onary models was com puted with the YREC7 in its 
non-rotating configuration (jGuenther et al ., 1992). All models are evolved from 
fully conv ective pre-main sequence ( PMS) to a stage of subgiant. The OPA L 
EOS200 1 dRogers fc Navfonovl l2002f) OPAL opacity (|lglesias k Rogerslll996l ). 
and the lAlexander k Ferguson! ( 19941 ) opacity for lo w temperature were used . 
These opacity tables have the solar mixtures given bv lGrevesse k Noels! ( 19931 ). 
The m odels took in t o acco unt diffusion of helium and metals, using the prescrip- 
tion of lThoul et al.l (| 19941 ) . Energy transfer by convection is treated according to 
the standard mixing-length theory, and the boundar ies of the convection zo nes 
are determined by the Schwarzschild criterion. See iDemarque et al.l (|2007| ) in 
details for the YREC. 



3.2. Computational method 

The position of a stellar model in the Hertzsprung-Russell (H-R) diagram 
depends on five modeling parameters: the stellar mass M, the mixing- length 
parameter a, the age of the star t, the initial hydrogen abundance Xi and 
metallicity Zi- For /i Her we take the ratio of the heavy-element abundance 
to hydrogen abundance on the stellar surface as an observable. Thus we have 
three observables (L, T e /f, (Z/X) s ) and five unknowns (M, a, t, Xi, Zi). 

In order to reproduce the observational constraints, we construct a grid of 
models with various masses, initial element abundances and mixing-length pa- 
rameters. For each stellar model, low- deg ree p-mode frequ encies are computed 
using Guenther adiabatic pulsation code (|Guentherj . Il994l ). To find the set of 
modeling parameters [M, a,Zi, Xj, t) that leads to the best ag reement with the 
observational constraints, following lEggenberger et al.l (|2004al ). we perform a \ 2 
minimization. The function \ 2 is defined as follows 



X 



Xclas 



x 2 ■ 

Aoscz ? 



(1) 
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where 



and 

W ,,theo 

,2 



(2) 

1 " ,,tneo obs / r) \ 



1=1 

Here the tr()'s are the errors on the corresponding observations, N is the number 
of observed frequencies, a — 1.8 //Hz is the resolution on the observed frequen- 
cies, and (D v ) — X)i=i( i 4' leo — v° hs )/N . The model which can minimize Xcias 
and Xosci at the same time will be considered to be the best one. 



4. Results 

Firstly, we computed a grid of evolutionary tracks for models with masses, 
mixing length parameters, initial metal and hydrogen mass fractions respec- 
tively in the ranges of 0.86 M© < M < 1.16 M©, 1.5 < a < 2.10, 0.025 
<Zi< 0.033, 0.57 < X, < 0.70 with a primary resolution SM = 0.02, 5a = 
0.2, 8Z — 0.002, and SX — 0.02. Computational results show that the evo- 
lutionary tracks for models with 0.9M Q < M < 1.15M Q and 1.5 < a < 2.10 
approximately span the error box in [T e ff, L, (Z '/ X) s ]. This is because the fact 
that a decrease of the mass can be compensated by a decrease in hydrogen and 



meta l abundances to get the same position in the H-R diagram (jMeng et al 



20081 ). Moreover, the evolutionary tracks imply that /i Her is in the post-main 



sequence phase of evolution. In this phase, the tracks are almost parallel to the 
T e //-axis (see Figs. [T][2]). Additionally, a variation of the mixing length param- 
eter mainly changes the radius , but has almost no influence on the luminosity 
(jKippenhahn fc Weigertl Il990h . Therefore, with increasing the mixing length 



parameter a, the evolutionary track moves almost horizontally to the left of H-R 
diagram and vice versa. Thus the models with various mixing length param- 
eters can span the same position in the H-R diagram (see Fig. [5]) at different 
ages. Table[3]gives the characteristics of three models, M2a, M2 and M2b, with 
different mixing length parameters but a same surface metallicity and position 
in the H-R diagram. These models have different xtsti' 1 - e - different oscillation 
frequencies, reflecting the differences in the internal structure of the models; 
for example, the differences in the central helium-core mass and density. This 
indicates that non-asteroseismic observational constraints do not enable us to 
determine the mixing-length parameter a for \i Her, an evolved solar-type star, 
but asteroseismic observations could provide a constraint on the mixing length 
parameter. 

For the sets of modeling parameters that lead to agreements with 
the observational constraints, we calculated the models with a fine 
resolution 5M = 0.01, 6a = 0.05, 6Z = 0.0005, and 5X = 0.005 in the 
vicinity of the parameters. We obtained many models that can almost 
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Figure 1: Evolutionary tracks in the H-R diagram for models of /i Her. Left panel: full tracks. 
Right panel: close-up of H-R diagram in the vicinity of /i Her. The positions in the H-R 
diagram of the models listed in Tableware indicated. The box shows the position of /i Her. 



Table 2: Models for /i Her. The superscripts or subscripts of some modeling parameters show 
confidence limits. 





logL/L© 




0.431 ± 0.025 




Observational 


Teff [K] 




5500 ± 90 




constraints 


{Z/X) s 
Av 




0.040 ± 0.006 
56.50 ± 0.07 








Ml 


M2 


M3 




M/Mq 


l.lO-o.oi 


1 00+ 01 


0.93 ± 0.015 


Modeling- 


a 


1.75 ± 0.25 


1.75 ± 0.25 


1.75 ± 0.25 


parameters 


Zi 


29+0.001 


0.0275t°;^ 


0.026 ± 0.002 




Xi 


0. 6 5o+o.oo5 


0.6051°;^ 


0.572 ± 0.01 




log L/Lq 


0.4550 


0.4386 


0.4307 




Teff [K] 


5454 


5494 


5492 


Model 


Z s 


0.0272 


0.0257 


0.0245 


characteristics 


(Z/X) s 


0.040 


0.040 


0.041 




R/Rq 


1.894 


1.831 


1.816 




age [Gyr] 

X.osci 


6.556 ± 0.03 


6.433 ± 0.04 


6.303 ± 0.03 




1.05 


1.07 


4.55 




Xdas 


1.18 


0.10 


0.03 




x 2 


2.23 


1.17 


4.58 






56.44 


56.58 


55.34 




Su 02 


5.00 


4.94 


5.09 




{D u ) [//Hz] 


28.10 


31.20 


-0.53 




He-core mass [Mq] 


0.100 


0.101 


0.105 
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Figure 2: Evolutionary tracks in the H-R diagram for models with Z,; = 0.0275, Xi = 0.605, 
and M = 1.00 Mq, with the lines differing in the value of the mixing-length parameter, a. 
The box shows the position of p, Her. Open circles indicate the locations of the M2a, M2 and 
M2b models listed in Table |3] 



Table 3: Characteristics of the models illustrated by open circles in Fig[2] 





M2a 


M2 


M2b 


M/M & 


1.00 


1.00 


1.00 


Zi 


0.0275 


0.0275 


0.0275 


Xi 


0.605 


0.605 


0.605 


a 


1.60 


1.75 


1.90 


age [Gyr] 


6.374 ± 0.04 


6.433 ± 0.04 


6.471 ± 0.03 


log L/Lq 


0.4387 


0.4386 


0.4387 


Teff [K] 


5494 


5494 


5493 


(Z/X) s 


0.040 


0.040 


0.041 


R/Rq 


1.832 


1.831 


1.832 


0.099 


0.097 


0.101 


x 2 

Aosa 


3.50 


1.07 


3.37 




56.65 


56.58 


56.46 


lo g/°c [g/cm 3 ] 


3.636 


3.800 


3.942 


He-core mass [Mq] 


0.091 


0.101 


0.109 
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minimize the Xosd an< ^ within the observational error box. Fig. [3] 
shows the xtsd as a function of mass and age. From this figure we can 
see that the models with mass 1 - 1.1 Mq and age 6.2 - 6.7 Gyr can 
better reproduce the observed frequencies. We also constructed the 
models with higher resolutions of modeling parameters than the fine 
resolution. But results are not sensitive to the resolutions. Performing 
the x 2 minimization described above, we found a solution: M = 1.00 Mq, a 
= 1.75, Z, = 0.0275, X, = 0.605, and t = 6.433 Gyr, marked M2. Tabled] 
lists characteristics of this model. The confidence limits of each modeling 
parameter correspond to the maximum/minimum values it can reach 
when other parameters are fixed, in order that the generated models 
fall within the observational error box. Corresponding evolutionary tracks 
are shown in Figure [TJ Although the value of Xcias an d Xosd °f tne M2 model 
is not the lowest one respectively, this model is almost able to minimize Xcias 
and Xosd a t the same time and has a lowest x 2 ■ Fig. [3] shows the model 
with a mass of 1.04 Mq has a lowest Xosd- But f ne va l ue of x 2 of this 
model is larger than that of M2. The mean large and small separation of 
the M2 model is 56.58 and 4.94 /^Hz, respectively, which are in good agreement 
with those obse rved. However, th e ma ss of M2 mode l . 1.00 Mq, is less than the 
values given by Fuhrmann (1998) and Takeda et al.1 ( 20051 ). They showed that 
the mass of /i Her is about 1.14 Mq. Table also gives a model Ml which has 
a mass as high as 1.10 Mq. This model is able to minimize Xosd m * ne a S e 
of 6.556 Gyr but has a higher value of Xcias- This implies that increasing the 
hydrogen abundance and mass at the same time for a fixed value of the mixing 
length parameter leads to an increase in the ~x 1i„„ in order to attain the sam e 
value of Xosd- This scenario was also found bv lEggenberger fc Carrier ( 2006h . 
The mean large and small separation of the Ml model is 56.44 and 5.00 /iHz, 
respectively, which are almost same as those of the M2 model. 

In Fig. 31 we plotted the differences between the observed and theoretical 
frequencies for the Ml and M2 model. Moreover, the theoretical and observed 
frequencies are compared by plotting echelle diagrams of the Ml and M2 mod- 
els in Fig. [5j The systematic difference {D y ) between computed and observed 
frequencies has been considered in these figures. The systematic difference {D v ) 
is 28.1 and 31.2 /iHz for the Ml and M2 model s, respectively. This system- 
atic difference also exists in oth er stellar models ( Christensen-Dalseaard et all 
19951 : lEggenberg er et all lioolal lbl: lEggenberger fc CarrierM2006h . Fig. |4] shows 



that the differences between the observed and theoretical frequencies are very 
similar for both models at high frequency, however they are different at lower 
frequency. Compared to the M2 model, the Ml model badly reproduces the 
observed frequencies of 1034.0 and 1061.2 /iHz, whereas the M2 model does not 
reproduce the observed frequencies of 947.6 and 1124.8 /iHz well. These com- 
parisons between individual observed and theoretical asteroseismic frequencies 
for the Ml and M2 models do not allow us to differentiate which model is better. 
Additionally, both models have the almost same values of the mean large and 
small separations reflecting a similarity of the internal structure between the Ml 
and M2 models. For example, both models have an almost same helium-core 
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1 1.1 

Mass (M ) 



6.5 7 
Age (Gyr) 



7.5 



Figure 3: Lowest values of X oac i as a function of mass and age. Open circles indicate models 
in accordance with the classical observational constraints, while dots show models that do not 
fall within the observational error box. 




1000 



1200 



1400 



v OitHz) 



1600 



Figure 4: Differences between computed and observed frequencies. Crosses and filled circles 
correspond to the Ml and M2 model, respectively. The systematic difference (D v ) has been 
taken into account in this figure. 
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1=2 t=0 
I 



1=3 t=\ 



20 40 60 

(f - <D„» modulo 56.5 (yuHz) 



1=2 1=0 
I 



1=3 l=\ 



20 40 60 

{v - <D„)) modulo 56.5 (/xHz) 



Figure 5: Echelle diagrams f or the Ml and M2 model, respectively. Filled triangles refer to the 
observed frequencies {jBonanno ct al., 2008), while open symbols correspond to the theoretical 
ones. Squares are used for modes with I = 0, triangles for I = 1, diamonds for I = 2 and 
circles for I = 3. 



mass. 

However, Fig. \5\ shows that both the Ml and M2 models badly reproduce 
the observed frequencies of I = 1. We searched the models which can better 
reproduce the observed frequencies of I = 1. A model labeled M2b, which 
was shown in Table El was found. However, the value of Xosd f° r the M2b 
model is larger than those of the Ml and M2 models, which is mainly due 
to differences in the modes of I = 1 of these models. If the contribution to 
xlsci from the I = 1 mode near 1238 /zHz is neglected, the value of Xosci f° r 
the M2b model will decrease from 3.37 to 0.72. Echelle diagram for the M2b 
model is shown in Fig. [6l Except for the I = 1 mode near 1238 /iHz, the M2b 
model reproduces the observed frequencies well. Moreover Fig. [5] shows that 
the observed and theoretical frequencies with I = 1 have a similar behavior: a 
zigzag echelle diagram. The modes of / = 1 show up a deviation from their 
expected asymptotic values. Fig. [7| also shows that the large separations for 
modes of I — 1 deviate from their expected asymptotic behaviour. This may 
be because that the formation of the central helium core leads to an increase in 
the frequencies of the g-modes. 'When the frequency of a g-mode approaches 
that of a p-mode, the tw o modes undergo an avoided crossing, where they 
excha nge physical nature' (jAizenman et all ll977tlChriste nsen-Dalseaard et al 



1995h . Figs. [6]and[7]show that the I = 1 modes of the M2b model between 1200 



and 1300 /iHz may undergo the avoided crossing. Noting the helium-core mass 
of the Ml and M2 models are almost same, but that of the M2 and M2b models 
are different, the behaviour of the I = 1 modes might be very sensitive to the 
internal structure of the stars. It may be applied to extract the information of 
the helium core. 

We also considered a Xosd which does not contain (£>„). For models falling 
within the observed constraints on effective temperature, luminosity and surface 
metallicity, we computed this new Xosci- A solution with M = 0.93 M , a = 
1.75, Z, = 0.026, Xi = 0.572, and t = 6.303 Gyr was found. The evolutionary 
track and the position of this model (marked M3) in the H-R diagram are also 
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(v - (D„» modulo 56.5 (jiHz) 



Figu re 6: Echelle diagrams for the M2b model. Filled triangles refer to the observed frequen- 
cies llBonanno et all 12008), while open symbols correspond to the theoretical ones. Squares 
are used for modes with I = 0, triangles for I = 1, diamonds for 1 = 2 and circles for / = 3. 
The systematic shift (D v ) is 28.4 /iHz for this model. 



shown in Fig. [TJ The characteristics of the M3 model are given in Tabled The 
systematic difference {D v ) between computed and observed frequencies is -0.53 
/xHz for this model. However the value of the new Xosd 1S 4-55, which is much 
larger than that of the Ml and M2 models. Moreover the echelle diagram of the 
M3 model is plotted in Fig. [8] This echelle diagram shows that the M3 model 
can reproduce the frequencies in the range of 1034 < v < 1437 /iHz but badly 
reproduces the modes at the lower and higher frequencies. Moreover its mean 
larger separation is only 55.49 /xHz. Thus this model is in disagreement with 
the asteroseismic constraints. 



5. Discussion and Conclusion 



The mass of 1.00 M p for the M2 and M2b models is less than the value in 
the previous literatures ( Fuhrmann L 1998; Tak eda et all 20051) . w hich may result 
from that the results of Fuhrmann ( 19981 ) and Takeda et all ( 2005 ) were obtained 
without the asteroseismic constraints. The initial hydrogen mass fraction is 
0.605 and 0.650 for the M2 and Ml model, respectively. Although there is a 
difference of 0.0015 in Z between the Ml and M2 model, the increase of the 
hydrogen or decrease of helium is mainly compensated by an increase of the 
mass in ord er to reach the s ame l ocation in the H-R diagram: a helium-mass 

" 1993f) . The Ml and M2 models have the same 

This implies that 



degeneracy ([Lebreton et al 
value of Xosd an d the almost same asteroseismic features 



the helium-mass degeneracy is difficult to be removed, even if the asteroseismic 
constraints are taken into account. 
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Figure 7: Large separations of theoretical frequencies of different models. The solid line shows 
the large separations of the 1 = modes, while the dashed line indicates the large separations 
of the I = 1 modes. 
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Figure 8: Same as Fig. [5] but for the M3 model. Left panel shows frequencies modulo the 
value of observed large separation, while right panel presents frequencies modulo the value of 
theoretical large separation. 
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Comparing with the observed frequencies, the theoretical frequencies have 
a systematic shift. The mode of I = 1 near 1238 fiRz deviate from its 
expected asymptotic value in Fig. [6j It may undergo an avoided 
crossing and be a mixed-mode. The mixed-mode could have higher 
inertia than other p-modes. The shift of this mode may be different 
from that of other p-modes. Thus we divide the Xosd into two parts: 
one Xoscm f° r mixed-modes and one Xoscp f° r & U modes except for 
mixed-modes. Assuming only the I = 1 mode near 1238 /iHz is a 
mixed-mode, we obtained the value of this new Xosd * s 1-07, 1.10, 
3.61 and 0.72 for Ml, M2, M2a and M2b, respectively. Under this 
assumption, M2b has a smallest Xosd- 

We confirmed the results that an analysis of the H-R diagram does not allow 
us to determine the m i xing l ength parameter for an evolved solar-type star 
( Fernandes & Monteirol . 20031 ) but the observed oscillation frequencies could 



provide a constraint on this parameter. The internal structure of the evolved 
solar-type stars is sensitive to the mixing length parameter at given (T e f f,L). 

In this work we constructed the models for the \i Her using the Yale stellar 
evolution code. By combining the non-asteroseismic constraints with the as- 
teroseismic observations, we find that a model for fi Her can reproduce the all 
non-asteroseismic and asteroseismic constraints well: the model with a mass 
of 1.00+° 02 M ®i an a S e * = 6.433 ± 0.04 Gyr, a mixing-length param- 
eter a — 1.75 ± 0.25 , an initial hydrogen abundance X* — 0.605+q'oq 5 
and metal abundance Zi — 0.0275+° ) '° ) ° ) 'j ! . However, the models with 
mass 1 - 1.1 Mq and age 6.2 — 6.7 Gyr also can reproduce the non- 
asteroseismic and asteroseismic constraints. Existing observational 
constraints do not rule out those models. 

The modes of I = 1 show up the characteristic avoided crossings, which 
may be applied to test the internal structure of an evolved solar-type star. 
The asteroseismic observations put important constraints on the models for 
[i Her, but they are not enough to really test the differences in the models. 
More accurate oscillation frequencies, especially the modes of I = 1, are need to 
investigate the internal structure of this type star. 
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